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Abstract 

The condition for = to be an eigenvalue of the operator \/— A + — m + XV 
is obtained through the use of the Birman-Schwinger principle. By setting E = —o? 
and using the analyticity of the corresponding Birman-Schwinger kernel, the series 
development of (A~^)(a) is obtained up to second order on a. 
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1 Introduction 



The Birman-Schwinger principle [Il[2] has been extensively used in several branches of physics 
and mathematics. The original motivation for its implementation was the possibility of 
having eigenvalue counting functions for the Schrodinger equation. In this approach, the 
potential term is taken to be of the form XV with A being a constant and V a potential 
satisfying a certain number of properties (e. g., is non-positive). The eigenvalue counting 
functions thus obtained will depend on A, and in general, for a given operator O, they can 
be defined by the following expression 

A^(A) = {number of fij G Spec(0),/ij < A}. 

Research on this class of functions was actually initiated by Weyl for the Laplace operator 
A [31 S], with an emphasis in the following years on the asymptotic expression for A^(A) 
when A — >^ oo. More recently, the use of Microlocal Analysis has indeed allowed a more 
detailed study of asymptotics with some proposals to replace eigenvalue counting functions 
by spectral shift functions [5]. 

Generalizations of the method itself have also been proposed to deal with certain physi- 
cal situations arising in quantum mechanics [6j and modified versions have been formulated 
for example in phonon physics where perturbations on acoustic and electromagnetic waves 
can be described through the presence of exponentially localized eigenmodes of their corre- 
sponding operators [7, 8J, and where an analogy with the appearance of bound states for the 
Schrodinger equation can also be established. 

In [9| the eigenvalue equation 

+ XV)^ = (1.1) 

was studied through the use of the Birman-Schwinger principle and some conditions express- 
ing the fact that the solutions to the Schrodinger equation are in £^(M^) were obtained. A 
related study was also done for the Dirac operator in [10]. In both cases the potential V 
were assumed to satisfy a certain set of properties. In this work we shall apply the same 
techniques to the case of the pseudo-differential operator a/— A -|- m? — m, which was studied 
by Herbst in the presence of a Coulomb potential [H]. This operator had also been used in 
other contexts such as hadronic physics [121 [El El- For a discussion of several aspects of 
this relativistic operator and its relation to the Dirac operator we refer to [TSl [TBI [T4l [TT] . 

In the next section we mention some basic facts regarding the Birman-Schwinger princi- 
ple. Section 3 deals with the analyticity of the Birman-Schwinger kernel for VA + m? — m 
and the series expansion of A~^ up to second order on a. The Fourier transform of the 
coefficient proportional to is then presented in Section 4. Several expressions used in 
the main text are derived in Appendix A and Appendix B. Appendix C contains a simple 
estimate on A + — E)~^{x)\. 

2 The Birman-Schwinger principle 

For convenience we recall the steps leading to the Birman-Schwinger principle and how it 
can be used in order to extract some information about eigenvalues. One begins from a 
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general eigenvalue equation of the form 

Hi; = {Ho + XV)iIj = Eip (2.1) 

with Hq an operator representing the kinetic energy. We shall suppose in the following that 
the potential V e C^(M^) and takes negative values. The above equation then can be recast 
in the form 

ij = XiHo-E)-'\V\^p. 
By multiplying this equation by \V\^^'^ and defining := one obtains 

XK(f) = (j), K = \V\^/\Ho - E)-^\V\^/\ (2.2) 

The physical meaning behind these rewriting is that the operator K has an eigenvalue Aq ^ if 
and only if H has an eigenvalue E = Eq for X = Xq, the multiplicity in both cases being the 
same. One can think of this in some sense as transformation which preserves the multiplicity 
among the eigenvalues of these two different operators. 

It is clear that the main feature of the method is to interchange the role played by 
the constant A and the energy E: the usual approach in quantum mechanics when dealing 
with an equation as Eq. (11. II) is to consider the potential term as being a perturbation and 
under suitable assumptions, the techniques of perturbation theory can then be employed to 
obtaine a (convergent) solution for the eigenfunctions and the eigenvalues E in terms of 
the perturbation parameter A. The Birman-Schwinger principle on the other hand allows 
one to consider the energy E as a free parameter on which a series expansion for A can be 
obtained. In fact, for most practical reasons it is useful to set E = —o? and consider a 
series development on powers of the parameter a (which is identified with the fine structure 
constant) as 

(A"^)(a)"^ = Ao^ + aa + 6a^ + ■ ■■ . (2.3) 

After this is done, the question is to know if one can solve the above relation for a in terms 
of A and if the resulting series obtained in this way is analytic. As outlined in |9j, most of 
the steps leading to Eq. (12. 3p pose no difficulties on this regard but a careful analysis has to 
be done for the second part. 



3 Kernel in coordinate space 



In this section we study some properties of the Birman-Schwinger kernel for the operator 
V— A + m} — m. The integral kernel is given by Eqs. ( lA.lll and ( 1A.2I) in Appendix A, namely 



K{x,y) = \V{x)\'l^ 



m 



4:7t\x — y\ 



'ti\x-y\ 



H — E{m\x — y\] jj) 



TT 



\V{yt/\ (3.1) 



The parameter /i is determined from 



(m + Ef 



-2mE — E^. By setting 



E = — a^, one has /i^ = 2ma^ — a^; this will be used at the end of this section and in the 
next one. 



In order to use perturbation theory on this operator, it is necessary to proof that it 
depends analytically on the parameter /i. For convenience we shall set m = 1 in the following. 
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First, if one considers the case of real /i, it is more or less clear that Eq. (13.ip makes only 
sense as long as < < 1, and thus one can try to look for analyticity inside the disc < 1 
in the complex plane (we shall see that this is not so an arbitrary choice for the problem). 
Let us then consider {f,Kg) for f,g E C? . There are 4 different integrals over x to 
be considered: 

• (1 - |y|i/2e-'^l-l /"' ' cosh{i2z)Ko{z)dzV\'^^g) 

Jo 

/oo 
e~'^^K,iz)dz\V\'/'g) 
, 1 



To deal with them we shall argue in the same way as in [TO] for the Dirac operator. The 
argument is based on writing (/, Kg) as a sum of integrals of the form 



c(/i) // f{x)-^^e-^\^~y\giy)d'xd'y 



with h{x, y) a function in £°°(R'^ xR'^) having compact support and c(/i) an analytic function 

\x-y\'- 



of IX. By using Sobolev's inequality one then finds that f{x)l^^g{y) e ^^(M^ x M?). The 



fact that K as an operator can be bounded by operators with compact support in means 
that one is dealing with a meaningful object. In this way the whole problem of an analytic 
extension on n is contained in the expression of the coefficients c{jji) and the exponential 

For the first integral one immediately has c(/i) = (1 — fx'^Y^'^ and h{x, y) = |V^(x)|^/^|a; — 
j/| In the second case obviously there is no dependence on n and it suffices to note 

that with the function h{x,y) = \V{x)\^^'^\x — y\Ki{\x — y\)\V{y)\^^'^ , which has compact 

3 XX TB>3\ v,^^ fl^\ H^,y) ^ rll'TO'S X TtPS^i 



support and is in £°°(]R^ x M^), one has fix) , % g{y) G 



The third and fourth integral require more care: first, for the third integral one has to 
consider c(/i) = (1 — /i^), which is analytical on /i, and a function 

h,{x,y) = \V{x)\^/^\x-y\ /' '\os\i{^iz)K^{z)dz\V{y)\^'\ 

Jo 

However, it is clear that h^{x,y) < Fi(jj,)\V{x)\^^'^\x — y\\V{y)\^^'^ where 

POO 

= / cosh.{^z)KQ{z)dz 
Jo 



as an elementary calculation shows. This also points to the necessity of considering ^ 1. 
One can now apply Sobolev inequality to the product — l/| obtaining 
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thus that f{x) l^i^^p G £^(M^ x R^). Furthermore, by writing cosh(yuz) in terms of 
exponentials functions, it is clear that h^{x,y) is analytical on /i. 

Consider now the fourth integral: one has as before c(/i) = (1 — /i^) and 

/■oo 

h^{x,y) = \V{x)\'/^\x-y\smh{f,\x-y\) / e-^^' Ko{z)dz\V{y)\^/\ 

One can now see that h^{x,y) < ^c\V{x)\^^'^\x — y\\V {y)\^^'^ , where c = KQ(z)dz = 7r/2, 
using the fact that e~'^^ < e~^'^ for z G [a, oo) and that KQ{z)dz < KQ{z)dz for a > 0. 

Sobolev's inequality then leads as in the previous case to f{x) ^^i^^p G £^(]R'^ x M'^). 
As before, the function is analytical on /x since it involves in a simple way analytical 
functions of fi in its definition. From the previous considerations it follows then that (/, Kg) 
is analytical on /i so that Theorem 3.12 in [H] can be used to deduce that K is analytic in 
the sense of norm. One should notice that the fact that one is considering a potential V 
with compact support is fundamental for this argument to work. 

Let us now proceed to finding under which conditions is an eigenvalue of Eq. ( 12.11) . On 
general grounds one can argue in the following way: given an operator O with inverse 0~^, 
consider a function solution to the equation 

lo(f) = /io0 (3.2) 

where Iq = \V\^^'^0~^\V\^^'^ . This equation is similar in structure to Eq. (12. 2p . Let us now 
define a function u through the relation 

u = 0-'\V\'/'(P. 

This function satisfies the integral equation 

u = 0-'\V\'/^(j) = fio'0-'\V\'/'\V\'/'0^'\V\'/^(j) = fi,'0-'\V\u, 

namely 

uix) = /xo^ / 0-\x,y)\Viy)\uiy)d'y. (3.3) 
^/supp V 

In the case of a potential V with compact support, let us say fi, one obtains a development 
in series for large of the previous expression as 

u{x) = fio'l^r' [ c,iy)\V{y)\u{y)d''y + fi^'\xr^ [ C2iy)\Viy)\uiy)d'y + . . . (3.4) 
Jn Jn 

and therefore 

/ c,{y)\Viy)\uiy)d'y = fio [ c,iy)\Viy)\'/'<P{y)d'y = (3.5) 

must be imposed in order to have u G £^(]R'^). In the above the functions ci{y) and C2{y) 
are determined from the series development of 0~^{x,y) in powers of l/\x\ for large 
We have assumed that the negative powers of appearing in this series development are 
integers, which is clearly seen from Eq. ()3.4p . This might be not necessarily the most general 
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scenario, but as long as the first term has a factor \x\~'^ with 7 < 3/2, as in the case we 
are considering, Eq. (13.51) will be obtained. It can be pointed out that by using the relation 
/Xq = one has in the sense of distributions the following equation 

{0-i2o^\V\)u={0 + Ho^V)u = (3.6) 

and thus the condition stated in Eq. (13. 5p is necessary and sufficient to guarantee that is 
an eigenvalue of the operator O + fi^^V. 

We now go back to our case of study, O = \/—A + — m. We have from Eq. (13.11) that 
when E = 0, 

u{x) = / \V{y)\u{y)d'y + ... (3.7) 



2vr/io \x\ jn 

for large |x| since ci{y) = 1 in Eq. ( 13. Sp . Therefore, according to the previous discussion, 
the constraint J \V{y)\^^'^(f){y)d^y = is a necessary and sufficient condition for to be 
an eigenvalue. One should notice that a similar condition exists for the Schrodinger op- 
erator. This can be understood on the grounds of the non-local behavior associated to 
^/—A + w? — m: non-locality is contained in the Bessel function Ki{z) appearing in the 
function F{m\x\]n = 0) in Eq. (1A.2I1 and for large distances (at least greater than m~^) it 



becomes highly suppressed. In this situation the main contribution turns out to come from 
a Schrodinger-like term [first term in Eq. ( 13.11) ]. 

One can also try to see what happens for small values of In this case there are two 
terms to lowest order to be considered, 

Ai= f ^AV{y)\u{y)d'y (3.8) 

Jn \y\ 

M = I hv{y)\u{y) r^^dzd^y. (3.9) 

Jn \y\ J\v\ ^ 



and 



Both terms are constants, independent of \x\. If supp V = Q contains the origin, then 
for a non-divergent value of u{\x\ = 0) to exist, these constants should have finite values. 
The constant Ai also appears when one considers a similar situation for the kernel of the 
Schrodinger operator. If the product \V{y)\u{y) behaves like ?/~^ with /? < 2 near then 
it is well defined (strict inequality means Ai = 0). On the other hand, A2 is more singular 
due to the presence of the integral involving the Bessel function Ki. The first integral (from 
\y\ to infinity) diverges as 1/y near the origin. This means that finite values are possible in 
this case when the product \V{y)\u{y) behaves like ?/~^ with /3 < 1 (again, strict inequality 
means a vanishing value for A2). 



4 Perturbation Theory and Fourier Transform 



In what follows we shall consider perturbation theory for Eq. (13.11) . From this kernel we can 
immediately write for a small the series development K = Lq + (2m) ^/^aA + 2ma^B + . . . 
where 



L,{x,y) = \V{x)\'l^ 



m 



A7t\x — y\ 



^ J m\x—y\ ^ 



dz 



\V{y)\ 



1/2 



(4.1) 
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A{x,y) = -^\V{x)nV{y)\'/', (4.2) 

and 

1 f 1 1 1 9 fm\^-y\ 

B{x, y) = |r(x)| 2 - ^1' - + - ^1' - 



2|a; — y| 



oo 1 pmlx—yl 



+ ^^-^ / iro(^)^^^^ + ^ / i^o(^)^'t^^ l^(l/)r/'- (4.3) 

As mentioned before, the presence of the second term inside the brackets in Eq. (14.11) is an 
indication of the non-locality of the Herbst operator. Furthermore, by going into Fourier 
space, one can see that the kernel of Lq is not positive definite, and for large values of k it can 
take arbitrary positive and negatives values. This is due to the strong divergence produced 
by the integral of the quotient Ki{z)/z, which when coupled with the factor — ?/| is of 
order 3. 

Let us now find an expression for A^^ = {(f), Kcj)) [see Eq. (12.21) ] where satisfies L^cf) = fi^cf) 
and II0II2 = 1. We have 

(A(a))-i =(0, Lo0) + (2m)^/2a(0, A(j)) + 2ma^{<p, Bcfy + . . . 

am^l^ f r ...... .V . (4-4) 



V2' 



TT 



The expression for the coefficient h is rather cumbersome as can be deduced from Eq. (|4.3I) 
and no definite statement on the sign of this coefficient can readily be established. But 
before with the analysis of it, let us review the general procedure used to write the energy 
as a function of the parameter A. Assuming /iq 7^ 0, one deduces from the expression 

{\{a))^^ = /io + aa + ha^ + . . . 

the relation 

Aq "^A(a) = 1 — Aoaa + Ao(Aoa^ — h)a^ + . . . , Aq := /ig "^5 

for a small. The next step is to invert this relation and to find a as a function of A; it is 
here where care should be taken. If a does not vanish then one has 



«(^) = -T^(^ - ^0) + ^^?rT^(A - Ao)^ + . . . 



and therefore 



i?(A) = -«(A)^ = --l-(A-Ao)^ + 



On the other hand, if a vanishes then 

Aq ^A(a) = 1 - Xoba^ + ca^ + . . . 
and it follows that there are two different expressions for a, 



a±(A) = ± 



^li-bV Ag(-6)3/2 
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One obtains then 



E±{\) = -«±(A)2 = - Ao) + . . . 



m 



Ki-b) 

and moreover, the series development is done in powers of (A — Aq)^^^. For O = V— A + m^ — 
the coefficient a is proportional to \V{y)\^^'^(j){y)d^y^ according to Eq. (14.41) and 

then it suffices to recall the discussion leading to Eq. (13. 7p to conclude that the (first) second 
series development corresponds to = (not) being an eigenvalue. For the first series, when 
a 7^ 0, it is possible to appeal to the implicit theorem function to deduce analyticity of a(A) 
at the point Aq and thus of -E'(A). When a = 0, analyticity does not hold in the Schrodinger 
case [9] because of an argument making use of the fact that — A+AV^ only admits non-positive 
eigenvalues for all A > (Theorem XIII. 11 in [19] applied to G C^(M^)). Unfortunately we 
have been not able to find a similar result for the Herbst operator in the three-dimensional 
case. 

We also remark that in the case when a = 0, if A approaches Aq from above by real 
values then b should necessarily be negative for a±(A) to be well defined. In the case of 
the Schrodinger equation it is known through an argument involving Fourier transform that 
b takes negative values when a vanishes [9], in consequence the energy approaches a value 
Eq = from below as A Aq . For the Herbst operator one can proceed in the same way, 
i.e., one can write 

B{x,y) = \Vix)\'/'B{x,y)\V{y)\'/' 

and therefore 

6 = 2m2(0,S0) = 2m2(0, |y|i/2i3|r|^/^) = 2m\f,Bf) = 2m^ j B{k)\f{k)\^d^k (4.5) 

with / = The sign of b is then related to the behavior of the Fourier transform B. 

We should also note that the fact that a = guarantees that the expression in momentum 
space on the right hand side of Eq. (14.51) is defined as — 0. 

As seen from Eq. (14.31) not all functions appearing in that expression have norm in £^(R^), 
and it is then necessary to use Fourier transform in the sense of distributions to find the 
Fourier transform B{k). The details of this are given in Appendix B, here we only write the 
final result of the calculations which reads as 



2m^ j B{k)\f{k)\'^d^k = 2m j {ni^B{ma))\f{ma)\^d^a 



2m 



2TTa^ An^a^ n 



6w'^ + 5w^ + 2 1 
+ 



:i + w2)5/2 o-2(l + ^i;2)l/2_ 



3 1 2w^ 

+ 



27rV3 (1 + ^2)5/2 27ra2 (1 + ^2)5/2 



|/(ma)|Va (4.6) 



where a = k/m^ a = \a\ E [0, oo) and w = 27rcr. Even though a positive term is present in this 
expression, its contribution is dominated by similar negative terms as one can easily verify 
from Eq. ( 14. 6p . Therefore one concludes that b is (conditionally) negative and that — > 0~ 
when A — ^ Aq . Moreover, it also follows from this expression that in the non-relativistic 
case {k = ma 0) the leading order contribution is given by —n^^a^^ — (27r'^)^V^'^, an 
analogous result to the Schrodinger case. 
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5 Conclusions 



It has been shown that having wave-functions of the operator ^/—A + — m in £^(]R^) 
is closely related to the vanishing of the coefficient a in the series development given by 
Eq. f l2.3p . This coefficient has the same form as in the Schodinger case as seen from Eq. (14.41) 
and its vanishing can be recast as a statement on the fact that = is an eigenvalue of 
the equation (a/— A + — m + \V)(f) = Ecj). Moreover, the coefficient b in Eq. (12.31) can 
be shown to consist of a Schrodinger-like part and more involved contributions due to the 
non-local nature of the operator under study. Fourier analysis however can be used to show 
that b is conditionally negative. The physical consequence of these facts is similar to the 
Schrodinger case, namely, the phenomenon of "coupling constant threshold". 

When considering the Schrodinger operator with a spherically symmetric potential V 
the above characterization can be explicitly formulated as a condition for the presence of 
s-waves as solutions of the equation (—A + V)"^ = and indeed, it is only the behavior 
of \l/ at infinity, where the potential V vanishes, which should be considered to obtain this 
result [9j. 

In the case of the Herbst operator a similar criteria should be also possible, however 
there is a simple feature that is present in one case and not in the other: the answer for 
the Schrodinger operator can readily be given because one knows precisely the expression 
for A in spherical coordinates meanwhile for the Herbst operator one should study instead 
an integral equation. The following auxiliary problem then arises: is it possible to write 
an appropriate integral equation for A\E' = from which one can recover the well known 
radial part r', r~'^^, / = 0, 1, . . . , of the spherically symmetric solutions? The answer is yes. 
The next natural step is to see if the same process can be applied to the integral equation 
obtained from (V— A + — m)\I' = 0. This is a subject that we hope to discuss in a future 
work. 
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Appendix A: Inverse Fourier transform 

In this section we would like to derive the expression for ( V— A + — E)^^ in coordi- 
nate space. Since one is lead to deal with tempered function, the Fourier transform in the 
generalized sense should be considered. Several ways exist to calculate the Fourier trans- 
form of tempered functions, such as the e-prescription, where it is customary to consider 
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after the calculations only those terms with no dependence on e in the limit e — 0. This 
method was widely used in the evaluation of propagators in the early days of QED. The 
Hankel transform [20] is other method introduced to deal with tempered functions. It is 
known that the e-prescription, the same being the case with the Hankel transform, is related 
to the notion of finite part, Pf. of a distributional pseudo-function [21] pp 20] acting on 
Schwartz's space S of infinitely differentiable and rapidly vanishing functions, together with 
its derivatives, at infinity. We shall also use the fact that for a distribution /, the relations 
< J^if), >=< /, > and < J^-\f), >=< /, > where (p e S hold [22]. 



Accordingly one can write 
1 



3 ^47^2^2 _|_ _ £J 



(j){p)dy =E 



--E 



4y^2p2 _|_ ^2 _ ^2 
1 



4y]-2p2 _|_ ^2 _ ^2 

+ (-A + 



(j){p)d^p + 



4y^2p2 _|_ ^2 _ J^2 



(j){p)(fp 



']r3 A/47r2p2 + m^[47r^p^ + iJ? 
Jr3 47r V + K 



■(l){p)(fp 



+ (-A + /i2 



1 



/4^2p2 _|_ j^2 [47;-2p2 _|_ ^2 
1 



-(j){p)d^p 



^4jl2p2 _j_ ^2^/^^2p2 _|_ ^2 



(t){p)d^p. 



where we have set /i^ = rn?—E'^. Therefore, from < J^{f), 4> >=< f, ^(0) >=< J-" ^{J-'{f)), >, 
one has in the sense of distributions acting on S, 



^ 4^2p2 _|_ ^2 _ 



X) = ET-^ 



1 



4y]-2p2 _|_ ^2 



[X] 



^4^2p2 _j_ 77j2^47]-2p2 _|_ ^2 



^4^2p2 _|_ ^2|^47^2p2 _|_ ^2j 



The first term is a standard example in textbooks on Fourier transforms or potential 
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theory and can be easily found: 
1 



/ 



(j){p)d p 



-(47r2p2+^2-)j 



dt(f){p)(fp 



p poo p 

Jm.^ Jo J«.3 



1 



(47rt)3/2 



(47r)3/2 y 



-tj.\x\ 



Atc\x 



■T{(p){x)d^x 



with (j) & S. In the above we have use Fubini's theorem and an integral representation of 
the Bessel functions K^i^z). For the second and third term one should calculate the inverse 
Fourier transform 

1 



-1 



^4^2p2 _|_ j^2['4^2p2 _|_ ^2j 

Using the fact that the usual properties of convolution of the (inverse) Fourier transform re- 
main valid for distributions [22], then the problem reduces to the convolution of JF^^ 



1 



with 



way as done before with the result 



x). This latter inverse Fourier transform can be found in the same 



(x) 



I 



^4^2p2 _|_ j^2 



, / X .•^ m f Ki(m\x\) ^, , , ,q 

(j){p)£p=— ^\ ^' T{(f)){x)d'^X. 



\x\ 



Therefore, in distributional sense, one has the following expression 



47^2^2 _|_ ^2 _ jt; 



E 



[X 



Airlxi 



-fi\x\ 



+ (-A + /i2 



Ki{m\ ■ I) 



47r| 



2772 



[X 



+ 



m[im? — jj^) \e '^''l Ki{m\ ■ |) 

— ; ;— * : : 



E 



Anix 
+ 



-^l\x\ 



[X 



+ 



m Ki{m\x\ 
2^ 



\x\ 



in-" ' ' 



[X 



It is clear from the above that this expression might be treated indeed as an ordinary function 
if the convolution in the last line makes sense. In fact, using the expansion [23] 



-XR 



XR 



-x:(2/+i) 

vr ^ — ^ 



TT 



2Ari 



^+i(Ari) 



TT 



2Ar< 



■^Hi(Ar2) 



P;(C0S7) 
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where i? = [r^ + — 2rp 005(7)]^/^, ri = min(r, p),r2 = max(r, p), only the term with / = 
is seem to give a non- vanishing contribution for the convolution after performing the integral 
over the angular variables. Therefore, one obtains 



e-^l-l Ki{m\ ■ I) 

— ; — ;— * : — : 



47r 2p 



TT 



TT Y 2p|a;| 2 
An 2p 



ni^)K,{y)ydy 



+ 



TT 



47r 



/i(p|a;| 



^K.{^^)My)ydy 



m?\x\ 



— smn — 

Hy \m 



+ 



4:71 sinh(p|x|) 
m p|x| 



Ki{y)ydy 
e-^yl"^K,{y)dy. 



The explicit expressions of the Bessel functions Ii/2{.x) and 7^^1/2(0;) have been used in order 
to write the above relation. The final outcome of the previous calculations is then 
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47r X 
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-^i\x\ 



H — F{m\x\; fi) 



IT 



1 - 



2 \ 1/2 



TT 



-F(m|x|; p) 



(A.l) 



with 



F(m|x|;p) =fri(m|a;|) 



— ;— * : : 



47r 



\X\ 



[X 



--Ki{m\x\)+ 1-^ 



p- 



rra|a;| 



sinh ( — ) Ki{y)ydy 



+ — sinh(p|x|) / e-^'y'"'Ki{y)dy 



(A.2) 



=iri(m|x|) + 1 - 



/ cosh ( — ] KQ{y)dy 
Jo \m / 



— sinh (p|x|) 



In the above we have used some relations among integrals involving Bessel and exponentials 
functions [24l [23] in order to write the last expression. 



Appendix B: Fourier transform of B{x, y) 



For convenience we present here the calculations leading to Eq. (14. 6p in the text. In the 
expression for B, several tempered functions in the sense of distributions appear and their 
Fourier transforms, which will be also tempered functions, can be found in this context. 
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In the following we shall apply the Hankel transform to a well-behaved function suited 
to the problem (strictly speaking the Hankel transform refers to the one-dimensional case) 
and then argue its validity for more general cases. 

To begin with we shall use the formula [20l §43] 

27r 



n-2 . 

k 2 Jq 



Here n is the dimension of the space considered. Let us apply this to the function f{x) 

-a rnikl 
10 



\x\ 



J™'"^' g{z)dz with < a, we have 
27r 



n — 2 

r 2 Jo 



drr 



dzg{zy/^Jn^{27Tkr) 



27r 



n — 2 

r 2 Jo 



271 



drr''^'^~''Jn^{2iTkr) / dzg{z) 
1 



1 1 



J 



dzg{z) / duu'^^'^ °Jn-2 (u) 

J z/m 



In the above we have defined a dimensionless parameter w = 2nk/m and assumed that the 
function g{z) is well behaved so that Fubini's theorem can be used. For our function g{z) 
we shall take indeed g{z) = z'^Kq{z),0 < [3. The previous expression is justified as long as 
n/2 < a < n since the integrals involving the Bessel function Jyiu) are then well defined. 
After using the relations 



/"OO 

/ s''J^{s)ds = 2^r( 
Jo 



and 



one arrives to 



•^0 fc=0 



2 ' ' 2 

(-l)'=(x/2)'^+2'= 



k\{fi + iy + 2k + l)r(z/ + + 1) 



2/3+r 



r(f; 



-w 



2-a 



i-w^r r(m + ^+"-°+^ )2' 

ml I'm, -^ Wilis' 



m=0 



m!(m+^) r(m+f) 
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The sum inside the brackets is just a hypergeometric function; the final result then reads as 



1 



2/3+n/ 



2-a-l 
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r(f; 



n-ao/^+w/^-g I 2 ^ 

(n-a)rm 



,n — a /3 + n — a + 1 /3 + n — a + 1 n n — a 
dA2' 2 2 2 2' 



-w 



(B.l) 



Using a reasoning due to L. Schwartz [25l pp 113], this expression can be seen to hold also for 
other values of a, in particular for a < 0, if one keeps in mind that the function / should be 
then understood as a tempered function as well as its Fourier transform. More care should 
be taken when n — a = —2h, a = —2h orP + n — a + l = —2h, h being a positive integer, 
where poles appear in the T function, but the appropriate modification to be used for those 
cases has also been pointed out by Schwartz. This modification will not concern us however 
because of the values of a {±1 or 0) appearing in Eq. (14. 3p and of the dimension n = 3. As 
a particular case then one can consider a = l,/3 = 0,n = 3, 



TT 



1/2 



W 



(27r)i/2 p 

1 1 
2P(l + w;2)i/2- 



2V2 2V 



^1(3/2)3^2(1, 3/2, 3/2; 3/2, 2; -w') 



This corresponds, taking into account a factor —2/{Tirn?) and replacement k = ma, to the 
second term inside the brackets in Eq. (14. 6p . 

It is also straightforward from the previous computation to identify the Fourier transform 



x\ 



of/(x) 



POO 

" / g{z)dz,0 < a, this is given by the second term in Eq. (IB. II ). namely 

J m\x\ 



' m\x\ 
1 



-W 



,n — a (3 + n — a + l (3 + n — a + l n 



X2F2K 



n — a 



-w 



(B.2) 



This formula can also be generalized to other values of a than those in the interval {n/2,n), 
care needed only to be taken in more detail at the points mentioned in a previous paragraph. 
Taking this into account, if we now calculate the case a = 0, /3 = 1, n = 3, the following 
expression is obtained 



1 



,5/2r(5/2)^ 



f27r)3/2 k3 



^^2^/2 ^,^^^^3^2 (3/2, 5/2, 5/2; 3/2, 5/2; -w^) 



31(3/2)' 

A ^^323/2 I 

4nk^ (1 + ^2)5/2 



corresponding, times a factor 2/(7rm) and replacement k = ma, to the positive term in 
Eq. diS). 
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Appendix C: An estimate 



We rewrite Eq. fjA.lll as 



(V-A + m2-E)-^( 



m 



Att\x\ 



( ^''^ 2 

\ ) IT 



We would like to show that for a fixed value of /i, each term inside the brackets is bounded 
as function of This is obviously true for /io(|a;|) = |x|e~'^l^l (< and also for the first 
term coming from the product |a;|F(m|x|; /i), namely = |x|i^'i(m|x|) (< 1). The next 



two terms to consider are 



h2i\x\) = Ixle-'^l^l cosh (^) Koiy)dy 

Jo ^mJ 



and 



/isda;!) = sinh (/i|a;|) 
h2 can be seen to be bounded as follows 



hilxl) < /i ^ / cosh (—) Ko{y)dy = 77- / cosh (—) [ exp -- - ^ 
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— ( 1 
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du , 
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2/i V 



^3/2 
2 \ -1/2 



Let us turn now to /is, one has 

1*00 1*00 
h^{\x\) < \x\e^\^\ / Ko{y)e-^y''^dy < \x\ / K^{y)dy 

J m\x\ J m\x\ 

< m\xo\'^ KQ{m\xo\) 

where m|xo| satisfies the transcendental equation KQ{y)dy = zKq{z) with solution z 
0.7451315. From this follows also that m|xo|-ft"o("^|a;o|) < 7r/2 and hence 

hsilxl) < 7t\xo\/2 < c7T/2m, c G [0.7451315, 00). 

Collecting all the previous results we have the following estimate 

_ [/ ,,2\l/2 

\{V-A + m^ - E)-Wx\)\ < 
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(C.l) 
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